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We consider a general model of branch competition that automatically leads to a critical branch- 
ing configuration. This model is inspired by the 4 — r\ expansion of the dielectric breakdown model 
(DBM), but the mechanism of arriving at the critical point may be of relevance to other branching 
systems as well, such as fractures. The exact solution of this model clarifies the direct renormalization 
procedure used for the DBM, and demonstrates nonperturbatively the existence of additional irrele- 
vant operators with complex scaling dimensions leading to discrete scale invariance. The anomalous 
" " | ", exponents are shown to depend upon the details of branch interaction; we contrast with the branched 

■ growth model (BGM) in which these exponents are universal to lowest order in 1 — v, and show 
that the BGM includes an inherent branch interaction different from that found in the DBM. We 

£N| . consider stationary and non-stationary regimes, corresponding to different growth geometries in the 

dielectric-breakdown model. 

Diffusion-limited aggregation (!]] produces complicated fractal structures by branch competition. Recently, by 
considering the dielectric breakdown model Q (DBM), a controlled renormalization group was developed within a 
4 — 77 expansion for this class of models ||. This expansion is based on considering an aggregate as a collection of 

, strictly one-dimensional branches; at rj = 4, the aggregate consists of a single branch. Given that the probability 
of a single branch pair surviving for time t is proportional to l/t u , to use the terminology of the branched growth 

\ model Q, a direct renormalization group was developed in an expansion in 1 — v = (4 — rf)/2. It was found that, 

■ without fine tuning, the system arrives at a critical point characterized by a scale-invariant tip-splitting rate so that 
the renormalized probability of a branch surviving for a time t is proportional to 1 jt and the rate of branch production 
and branch death balance at all scales. However, branching structures are common in other physical systems, and a 
similar 1/t kernel has been found in fracture systems ||. 

Thus, we will consider a more general model of branch competition, inspired by the 4 — r\ renormalization group. 
This model will be exactly solvable, clarifying the direct renormalization procedure employed for the DBM. We will 
find that the exponents depend on the bare kernel and on the details of branch interaction, becoming trivial as the 
v — > 1. Recently ||, it was shown that within the branched growth model (BGM) the fractal dimension depends 
only on v for small v\ here, this is shown not to be true forgeneral models of branch interaction, though the fact 
that the numerical values of the exponents within the 4 — r\ || and branched growth [^) expansions are similar is an 
indication that the BGM is a useful approximation. However, we show that the physics of the BGM involves certain 
assumptions which do not hold for the DBM, so that the BGM result is not exact for the DBM. The exact solution 
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of the general model will also reveal the presence of additional, irrelevant operators with complex scaling dimension, 

I 



a 

\ indicating the presence of discrete scale invariance, which has been argued to exist in DLA 

Finally, we will consider this model within two separate regimes. One regime will be analogous to the steady state 
regime in the cylindrical geometry in the DBM, while the other will be a non-stationary state, which corresponds to an 
initial condition of a single infinitely long vertical branch in cylindrical geometry, over time scales much shorter than the 
time required for the aggregate to reach the scale of the cylinder. We will find that the exponents characterizing local 
fractal dimension and mass-radius scaling are the same. This indicates that the difference between these dimensions 
observed in radial geometry j|] , as well as the non- trivial affine exponents observed in the early stages of growth in 
cylindrical geometry ]9| from initial conditions consisting of a horizontal line, are due not just to the non-stationary 
nature of the growth but also to the different geometry and initial conditions. Another motivation to consider different 
geometries is to compare the branched growth model computation, performed in a non-stationary regime, to the 4 — 77 
RG computation, performed in a stationary regime. 

A Model of Branch Competition — The dynamical state of the system at a given time t will be defined by a binary 
tree, with a set of times ti, one for each branch point in the tree, defining the time at which that branch pair was 
produced. There are two dynamical processes, tip-splitting which leads to the production of additional branches and 
branch competition. The first process is accounted for by assigning a rate g at which each tip splits, changing the 
topology of the tree by adding an additional branch point at that tip with time ti = t. See Fig. 1. 

The second process is accounted for by assigning a probability of a branch point i being removed from the tree 
due to competition of branches. When the branch point is created, one of the two branches below the branch point 
is randomly designated as the weaker branch, and the other as the stronger. When the branch point is removed, the 
weaker branch is removed from the tree. See Fig. 2. We pick the probability of removing branch point i to be 

v6{t-ti -1)— (1) 
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plus additional corrections due to branch interaction. To include the effects of branch interaction, we pick a phe- 
nomenological model for the interaction. We pick a scale factor x < 1, and declare that if on the weaker branch there 
exists another branch point with branching time exactly equal to t — x(t — tj) then the weaker branch will be removed. 
In this case, the lower branch point will have lived for a time x shorter than the branch point i. 

This model has the essential features of the 4 — r\ renormalization group. One difference is that within that procedure 
the details of branch competition had to be determined numerically by integrating the trajectories of several competing 
branches, where here they may be determined phenomenologically by the scale factor x. Another difference has to do 
with the time scale. Within the DBM, each tip has a given growth measure and the tip-splitting rate and probability 
of removing branch points are proportional to the growth measure at the tip normalized by the total growth measure 
for the cluster. Here, we instead take all tips to have the same growth measure (so that all tips have the same 
probability of splitting) and we do not normalize the tip-splitting rate. The lack of normalization simply changes 
the overall time scale; the choice of the same growth measure for all tips is simply taken to make the model more 
tractable analytically and does not alter the essential physics. 

In addition, in the 4 — r] renormalization group, branch points may also be removed as a result of interactions with 
parent branches of similar scale, rather than just as a result of daughter branches. We will consider this possibility 
later. 

Survival Probability — Let a branch pair be created at time t — 0. Define s(t) to be the probability that the branch 
pair survives until time t, assuming that the branch pair is not destroyed by removal of one of its parents. We find 

dts{t) = _ ^(t-l)s(t) g{1 (2) 

as g(l — x)s(xt)dt is the probability of creating a branch point at time (1 — x)t which survives until time t. The 
growth rules have been chosen so that Eq. (||) is exact; the process of creating and then removing a branch point does 
not effect the distribution of branches on the remaining stronger branch below that branch point. 

Searching for a scaling solution, suppose gs(t) = A/t. We find 1 — v = (1 — x)A/x. The constant A sets the 
probability, in time t, of producing a branch which remains in the tree for time i; this is the dimensionless (scale 
invariant) tip-splitting rate j|. As v —*> 1, A — ► 0. 

This result can also be obtained by a direct renormalization procedure || in which one expands the survival 
probability in powers of 1 — v and g: to zeroth order in g, gs(t) = A/t + A(l — v) logt/t + 0(1 — v) 2 . Then, to first 
order in g we find gs(t) = A/t + A(l — v) logi/t — A 2 (l — x) logt/ (xt), so that a fixed point of ts is reached only for 
the given value of A. 

To investigate the approach to the scaling solution, suppose instead gs(t) — j(l + /(f)). Linearizing Eq. (Q) about 
/ = 0, wc find 

df(t) 
d log t 



-(l-u)f(xt). (3) 



Eq. (|b|) is translationally invariant in r = logt, and has solutions / = e fer , where the eigenvalue fc is the scaling 
dimension. We find 

k = -{l-v)x k . (4) 

When v w 1, the eigenvalue with largest real part is fc = —(1 — v) + 0(1 — v) 2 . Then, for v w 1, where this 
eigenvalue is small we can make an approximation that f(xt) — f(t) and approximate the nonlinear problem by 
dts(t) — —vs(t)/t — g(l — x)s 2 (t)/x, which can be solved exactly for s(t),t > 1 as s(t) = (1 — v)/(g(l — x)t/x + ct v ), 
where c is an arbitrary constant. 

For l — i/ small, Eq. (|4|) has two solutions for real fc, as well as an infinity of solutions with complex fc. All these k 
have negative real part and describe irrelevant perturbations. For x 1 ^" = e~ x / e the two real solutions merge, and for 
x x ~ v < e _1 / e all solutions have complex fc. The presence of complex eigenvalues indicates that s(t) has an oscillatory 
behavior and that there is a discrete scale invariance in the corrections to scaling. The scale of this discrete scale 
invariance is in general not equal to x, so that it is not simply an artifact of the particularly simple form of branch 
interaction chosen, but rather a result of the fact that branches separated by a finite range of scales interact. As 
x x ~ v is decreased, one finds a fc with vanishing real part when fc = in/(2 logs), so x 1 ^ u = e~^l 2 . For x x ~ v < e _7r / 2 , 
there are complex eigenvalues with positive real part which describe relevant perturbations. In this case, we have 
found numerically that f(t) is driven to a new fixed point with non-decaying log-periodic oscillations, with scale not 
in general equal to x. The model considered here is useful for analyzing these effects, which are beyond perturbation 
theory in 1 — v , 
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Consider correlations: the probability, given that there is a branch point at time t\ remaining in the tree until time 
t'i, that there is another branch point at time ti remaining in the tree until time t' 2 . To study these, generalize the 
survival probability to a function s(t,t'), the probability that a branch point created at time t survives until time 
t + t'. We obtain the equation 

d t >s(t, t') = _ ini_zi>(vO _ g(1 _ X ) 8 ^ t / )s((1 _ x)t + xt ^ xt ,y (5) 

Defining gs(t,t') = y(l + f(t,t')) and linearizing we find 

flf /(*,0 = ~^Af(t + (1 - (6) 

This equation is translationally invariant in t and so we look for solutions f(t,t') = e M g(t'). We find dt>g(t') — 
— ^Z"' g(xt')e zlxt . For I — 0, this is the same as Eq. (||). Generally, for t' « Z _1 , we find the same discrete scale 
invariant solution as for Eq. (Q) as above. For i' >> the perturbations f(t,t') decay more rapidly due to the 
oscillations of the exponential. Thus, the two time scales t, t' are related, and we find discrete scale invariance in both. 

Different Geometries — We now consider the local "fractal dimension" of the cluster. The length of a branch is 
defined to be t — U, where U is the time at which the pair to which that branch belongs was created. The largest such 
t — ti sets a length scale for the cluster. The mass of a branch is defined to be the length of that branch, plus the sum, 
over all branch points j which lie within that branch at time t, of t — tj, plus the sum, over all branch points which 
were added to that branch and later removed, of the elapsed time between the addition and removal of the branch 
point. 

Define m{t) to be the average mass of a branch which is removed from the tree after a time t has elapsed since 
its creation (hence, this is a weaker branch). By averaging over the weights of sidebranches, we obtain the exact 
recursion relation 

t xt' 

m(t)=t + g J dt' J &t"m{t")(-d t »s{t")). (7) 
o 1 

Assuming a power law m(t) — t D , one finds a(a — 1) = Ax a ~ x . For v w 1, we find 

D = l + {l-v)-^ + 0{l-v) 2 . (8) 
1 — x 

Consider two different geometries. If we start the tree with a single branch and follow the dynamics above, this is 
analogous to starting the DBM with a single branch as a seed configuration and letting the cluster grow. In this case, 
Eq. (|8|) provides a scaling of the mass with the time. In another geometry, analogous to the cylindrical geometry, we 
modify the dynamics to always remove a weaker branch if t — ti > T, for some T setting a scale. In this case, after 
an initial non-stationary regime lasting for a time of order T, the mass of the cluster increases linearly with t, with a 
rate proportional to T D ~^: the mass of the largest branches, times the probability of producing such a branch. Thus, 
within this model the mass-radius scaling and local fractal dimension are the same up to a trivial difference of unity. 
We can generalize the model by including a possibility of removing a branch due to the presence of a parent branch 
of comparable size. Then, the probability of the earliest branch point i surviving till time t will scale as (t — ti)~ a ti h 
where a + b = 1. One will again find that the average time required to produce a branch point surviving for time t is 
of order t and the mass-radius scaling and local fractal dimension will again be equivalent. 

Comparison to Branched Growth Model — It has been shown Q that there exists a 1 — v expansion for the BGM 
similar to the 4 — r\ expansion for the DBM. One elegant feature of this expansion is that the lowest order fractal 
dimension is obtained without considering interaction of branches, but simply from the bare constant v. This seems 
surprising, as we have found within the model above, and within the 4 — rj expansion for the DBM, that unless we 
include branch interaction, the tip-splitting rate grows at large scales for v < 1, and a scale invariant fixed point is 
not reached. 

The resolution of this is that the BGM is defined in a way which inherently includes effects of branch competition. 
In the BGM, a branch is assumed to have a probability \/m u of surviving until it reaches a mass m, while within 
the model above the probability is defined in terms of the probability to survive for a time t. Now, given that a scale 
invariant fixed point can only be reached if a branch has probability 1/t of surviving for time t, then we must have 
the relation that 1/m" = 1/t, so that m = t 1 '", and the fractal dimension is 1/v. 

More formally, m(t) = t + Atlogt to order (1 - v) Q ,A l . We do not have an exact Eq. (^) for the BGM, but 
this equation is still correct to lowest order. One may still define a survival probability s(t) for a branch within the 
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BGM, and d t s(t) = -vd t m(t)s{t)/m(t) so s(t) = 1/t + (1 - v) log t ft - Alogt/t to order (1 - v) 1 , A 1 , and we find 
A = \- v = D - \. 

However, this particular form of branch competition in the BGM is not that found in the DBM, as found by 
numerically following the evolution of three branches. In some circumstances production of a daughter branch can 
actually reduce the competition of the parent branch with its sister || ; it is only when summing over all configurations 
that the increase in branch competition is obtained. Further, the competition of branches which is inherent in the 
BGM involves only an increased competition of a branch pair due to daughter branches. However, within the 4 — 77 
expansion, it was necessary to consider the interaction of branches with daughter and parent branches (which is in 
fact the strongest interaction numerically) to obtain the correct result. Thus, both the present toy model and the 
BGM are approximations to the physics near 77 = 4, although the BGM dynamics serves also as a good approximation 
for 77= 1 |§. 

Conclusion — We have examined a simple, solvable model for branching, finding a fixed point with a scale invariant 
tip splitting rate. The direct RG for the model is exact at lowest order. The DBM has similar behavior and a similar 
perturbative RG. The exact solution enables us nonperturbatively to find additional irrelevant operators leading to a 
discrete scale invariance. We contrast the behavior of this model, and the DBM, with that of the BGM, for which a 
similar 1 — v expansion is available. 
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FIG. 2. Removal of a branch below the circled branch point. 
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